EVALUATION  OF  A  PROPOSED  MODIFIED 
LOG-GAMMA  CONFIDENCE  BOUND  METHOD 
FOR  FLEET  MISSILE  SYSTEM  RELIABILITY 


Peter  Allen  Craig 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

EVALUATION  OF  A  PROPOSED  MODIFIED 

LOG-GAMMA  CONFIDENCE  BOUND  METHOD 

FOR 

FLEET  MISSILE  SYSTEM  RELIABILITY 

by 

Peter  Allen  Craig 

September  1979 

Thesis 

Advisor:                   M.  Woods 

Approved  for  public  release;  distribution  unlimited 


T18958 


UNCLASSIFIED 


SECURITY   CLASSIFICATION   OF   THIS  PAGE   rwhmr<  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


"\       «»OiT    MUlTiTS 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.   GOVT    ACCCUION    MO 


1.     RECIPIENT'S  CATALOG   NUMIO 


4.     TITLE  (end  Subtitle) 

Evaluation  of  a  Proposed  Modified  Log- 
Gamma  Confidence  Bound  Method  for  Fleet 
Missile  System  Reliability 


5.     TYRE  OF   REPORT   a   RERIOO  COVERED 


Master's  Thesis; 
September  1979 


4.     PERFORMING  ORG.    REPORT   NUMBER 


7.     AUTHOftfa.) 

Peter  Allen  Craig 


i.  CONTRACT  OR  GRANT  NbMlCMUj 


»   "ERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California   93940 


10.     PROGRAM   ELEMENT.  PROJECT     TASK 
AREA  *    WORK   UNIT  NUMBERS 


II.     CONTROLLING  OFFICE  NAME   ANO  AOORESS 

Naval  Postgraduate  School 
Monterey,  California   93940 


12.     REPORT   DATE 

September    1979 


IS.     NUMBER  OF   PAGES 


46 


14.     MONITORING   AGENCY   NAME   4    AOORESS//f  different  Irom  Controlling:  Olllce) 


IS.     SECURITY  CLASS,  (at  rftla  riport) 

Unclassified 


l»a.     DECLASSIFICATION/ DOWNGRADING 


SCHEDULE 


16.     DISTRIBUTION  STATEMENT  (<ti  thlt  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.     DISTRIBUTION  STATEMENT  (of  ine  mbmtrm<t  entered  In  Block  30.  II  dltleront  fm  Report) 


IB.     SUPPLEMENTARY  NOTES 


IS      KEY  WOROS  (Continue  an  reweree  tide  II  nocooemrr  end  identity  by  block  number) 

Fleet  Missile  System  Reliability 
Log-Gamma  Confidence  Bound 


20.     ABSTRACT  ( Continue  on  ravaraa  oldo  II  nocoocmry  and  Identity  by  block  number) 

A  statistical  method  is  evaluated  to  determine  its  accuracy 
for  estimating  lower  confidence  bounds  on  system  reliability 
of  a  mixture  of  missile  configurations  using  component  data. 


DO  ,  ^  1473 

(Page  1) 


EDITION  OF   <  NOV  ••  IS  OBSOLETE 

S/N    0  102-014- 6601   | 


UNCLASSIFIED 


SECURITY   CLASSIFICATION  OF   THIS  PAOE  (Wmm  Dele  tutored) 


UNCLASSIFIED 


«CUWT»   ClMtl»IC*tiow   Q0    Twit  »40»fW ■«..   n...   *»< 


#20  -  ABSTRACT  -  CONTINUED 

Monte  Carlo  simulations  are  used  to  establish  the  accuracy 

of  these  bounds. 


DD      ForrQ       1473 
S/N    (H&2-014-6601 


UNCLASSIFIED 


SKCU«iTV   eLAUl^lCATlOM   Or    Twif    •»«'»»►•   0«f«    !-.'.« 


Approved   for   public   release;    distribution   unlimited 

EVALUATION  OF  A  PROPOSED  MODIFIED  LOG-GAMMA  CONFIDENCE 
BOUND  METHOD  FOR  FLEET  MISSILE  SYSTEM  RELIABILITY 


by 


Peter  Allen  Craig 
Lieutenant,  U.  S.  Navy 
B.S.,  Marquette  University,  1973 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 
from  the 


NAVAL  POSTGRADUATE  SCHOOL 
September  1979 


ABSTRACT 

A  statistical  method  is  evaluated  to  determine  its 
accuracy  for  estimating  lower  confidence  bounds  on  system 
reliability  of  a  mixture  of  missile  configurations  using 
component  data.   Monte  Carlo  simulations  are  used  to 
establish  the  accuracy  of  these  bounds. 
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I.   INTRODUCTION 

A  statistical  method  has  been  proposed  which  obtains 
a  lower  confidence  bound  on  system  reliability.   It  is  a 
modified  log-gamma  procedure  developed  to  measure  fleet 
missile  system  reliability.   Monte  Carlo  simulations  were 
performed  to  evaluate  its  accuracy  as  an  estimate  for 
system  reliability.   Five  hundred  simulations  were  run  for 
each  of  twelve  cases  examined  at  80%  and  90%  confidence 
levels.   The  results  of  these  simulations  are  included  in 
this  paper.   Additional  simulations  were  performed  with  minor 
modifications  to  the  proposed  log-gamma  method.   These  changes 
are  documented  and  the  results  are  included.   A  comparison 
was  made  between  the  two  versions  on  their  accuracy  for 
estimating  the  lower  confidence  bound  on  system  reliability. 

The  reliability  equations  were  applied  to  a  hypothetical 
fleet  missile  system  configuration  and  analyzed  for  changes 
in  test  sample  sizes,   component  reliabilities  and  weighting 
factors.   The  proposed  procedure  was  determined  to  be 
significantly  inaccurate  for  small  and  large  amounts  of 
accumulated  test  data  on  missile  components.   It  also  has 
the  distracting  defect  that  larger  lower  confidence  bounds 
are  obtained  from  data  sets  with  one  failure  than  those 
obtained  from  data  sets  with  zero  failures. 


II.   MODIFIED  LOG-GAMMA  METHOD 

The  log-gamma  method,  in  its  more  general  form,  can 
apply  to  nonseries  as  well  as  to  fleet-mixture  populations. 
The  underlying  theory  is  contained  in  [Ref.  1] .   Examples 
of  cases  where  it  is  suspect  have  been  included  in  the 
following  chapter.   The  procedure  below  describes  the 
proposed  modified  log-gamma  method  as  it  is  applied  to  a 
series  system. 

Assume  that  in  a  series  system  there  are  k  components 
each  with  a  sample  size  n.,  where  i  =  1,  2,  ...,  k.   Let 
the  number  of  failures  be  f.  for  i  =  1,  2,  ...,  k.   Consider 
first  the  case  when  there  is  at  least  one  failure  in  the 
system.   Thus  Zf.  >  0. 

Let 


f  . 

R.   =   1 -  (2.1) 

in. 

l 


be  the  point  estimates  of  the  i-th  component  reliability. 
Then  the  equation 


k  „ 
R   =    tt  R.  (2.2 

i-1  X 


is  the  point  estimate  of  the  system  reliability.   Define 


R  =   R1/k  (2.3) 


and 

V   =   (1  -  R)    I      ~-  (2.4) 


1=1   l 


V  is  used  as  an  estimate  of  the  variance  of  -In  R.   It  is 
assumed  that  the  distribution  of  -In  R  can  be  approximated 
by  a  gamma  distribution  as  follows 


L-l   Lz/ln  R 
f(z)   =    — ,   z>0         (2.5) 

(:^-^)  r(L) 

-In  R 
where  z  =  -  In  R  and  L  and  ( — = )  are  parameters. 


Let 

*,  2 

(2.6) 


*      (-In  R) 


V 

and 

L   =   L*  +  2.25  (2.7) 

L*  is  the  method-of-moments  estimate  of  the  shape  parameter 
A  constant  term  2.25  is  added  to  L*  ,  the  shape  parameter 
estimate  in  the  proposed  modified  log-gamma  procedure.   The 
lower  (1  -  a)  confidence  bound,  R(l  -  a)  is  given  by  solving 
the  equation 

J^/X22l   a)  (2.8) 

R(l-a)   =   R      Zh,a 

2~ 
where  Xt£    ^s  the  lower  a-quantity  of  the  chi-square 

distribution  with  2L  degrees  of  freedom.   Interpolation  is 

required  if  2L  is  noninteger. 


If  there  are  zero  failures  in  the  system  (J  f.  =  0) 
let 


N*   =  -jf-  (2.9) 


i-1  ni 


where  N*  is  defined  to  be  the  effective  sample  size.   Then 
the  lower  1-a  confidence  bound  R(l-a)  is  computed  according 
to  a  binomial  confidence  bound  based  on  zero  failures  out  of 

N*7 — 

N*  trials  (i.e.,  R(l  -  a)      =    /a  ) .   If  N*  is  noninteger 
then  linear  interpolation  is  recommended  in  the  proposed 
procedure  but  it  is  not  necessary  because  the  same  formula 
could  be  used  for  N*  an  integer. 

The  modified  log-gamma  method  has  been  described  here 
for  both  zero  failures  and  one  or  more  failures  in  series. 
The  more  general  form  of  this  method  was  applied  to  an  actual 
missile  system  configuration  to  determine  the  lower  confi- 
dence bounds.   The  program  used  to  evaluate  its  accuracy 
has  been  included  in  Appendix  B.   The  complete  listing  and 
definitions  of  the  variables  used  in  the  program  are  listed 
in  Appendix  A.   A  description  of  the  more  generalized  method 
is  described  as  it  was  applied  to  the  specific  missile 
system  simulated. 

In  the  fleet  missile  system  examined  there  were  different 
groups  of  missiles  with  different  configurations.   The 
population  was  therefore  not  homogeneous  and  weights  were 
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assigned  to  the  different  groups.   There  were  14  components 
in  the  system  modeled  and  eight  mixture  weights  for  the 
subgroups.   The  input  data  consisted  of  f .  (the  number  of 
failures  in  the  i-th  component)/  n.  (the  sample  size  for 
the  i-th  component),  M.  (the  exponent  of  each  component) 
and  C.  (the  weights  applied  to  each  subgroup).   Point 
estimates  for  this  system  were  defined  as  follows 


(2.10) 


R. 

l 

= 

1       " 

f . 

i 

n . 

l 

Pr 

= 

5 

n 
i=i 

M. 

R.     X 

l 

/\ 

10 

M. 

PN 

= 

IT 
i=6 

R.     1 

l 

(2.11) 


(2.12) 


with  the  subgroup  reliability  point  estimates  being 


R     -   ^  Rll  R12 


/\  /\  /\ 


R      "   PN  Rll  R12 


/V  /N 


"  (3) 

R      -   PR  R13  R14 


/N  /N 


~(4)   -   pN  R13  R14  (2.13) 


/\     /\ 


~  (5) 

R      -   PR  R13  R12 


/\  S\ 


R    -  PN  R13  R12 
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•\  s\  /\ 


(7) 
R  -      PR  RU   R14 


/\  /\ 


R  ~      PN    Rll    R14 


and 


A  Q  A    It   i 

R      =         T       c.    R°;  (2.14) 

j-1       D 


The  variance  of  -In  R  is  then  estimated  by  V  given  by 
equation  (2.15) 


V   =   i    7    F   c.c.  R(l)  R(II)  S.  .        (2.15) 
«2   ..*■,_*•,   a.  I  li 


R    1=1  j=l 


where  S..  estimates  the  cov(z    ,  z  -*  )  and  where 
lj 

z     =   -  In  R    .   The  estimates  S.   are  found  by  solving 
the  following  equations. 


z.   =   -  In  R.  (2.16) 


14  14 

R   =   exp(-   [   M.z.  /   J   M.)  (2.17) 

.  i   li    .,   i 


i=l  *  *   i=l 


and 


VR 

= 

(1-  R) 

5      M.2 

1         L 
i-1     ni 

VN 

= 

(1-  R) 

10     M.2 

1         * 
.  L  r     n  . 

(2.1: 


(2.19 


i=6   i 


V.   =   (l-R)/n.     ,   i  =  11,. ..,14        (2.20) 
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Then  the  S. ,'s  are  solved  by  the  equations  listed  in  the 
program  in  Appendix  B  and  repeated  below. 


S(l,l) 
S(2,2) 
S(3,3) 
S(4,4) 
S(5,5) 
S(6,6) 
S(7,7) 
S(8,8) 
S(l,2) 
S(l,3) 
S(l,4) 
S(l,5) 
S(l,6) 
S(l,7) 
S(l,8) 
S(2,3) 
S(2,4) 
S(2,5) 


VR  +  Vll  +  V12 


=  VN  + 


Vll  +  V12 


VR  +  V13  +  V14 
VN  +  V13  +  V14 
VR  +  V13  +  V12 
VN  +  V13  +  V12 
VR  +  Vll  +  V14 
VN  +  Vll  +  V14 
Vll  +  V12 


=  V, 


=  0 


VR  +  V12 


=  V 


12 

V„  +  V 


11 


=  V 

=  0 

=  v. 

=  V 


11 


N 


S(2,6)  =  VN  +  V12 


S(2,7)  =  V 


11 


8(2,8)  =  VN  +  Vxl 

S(3,4)  =  V13  +  V14 

S(3,5)  =  VR  +  V13 

S(3,6)  =  V13 

8(3,7)  =  VR  +  V14 


S(3,8)  =  V 
3(4,5)  =  V 


14 


13 


S(4,6)  =  VN  +  V13 


S(4,7)  =  V 


14 


12 


S(4,8)  =  VN  +  V14 

S(5,6)  =  V13  +  V12 

S(5,7)  =  VR 

S  (5,8)  =  0 

S (6,7)  =  0 

S(6,8)  =  VN 

S(7,8)  =  Vu  +  V14 


(2.21) 


Finally, 


£   =   (-  ln  R>   +  2.25 
V 


(2.22) 
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and  DF,  the  degrees  of  freedom,  is  equal  to 


Thus 


DF  =  2L  (2.23) 


R(l-a)   =   R(DF/xDF,a}  (2.24) 
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III.   EVALUATION  PROCEDURE 


The  equation  for  system  reliability  is 


L      k    M. 

R    =    I   w.   n   p.  1  (3.1) 


j=l  J  i=l 


where 


L  =  number  of  subsystems 

w.  =  the  weighting  factor  of  the  j-th  subsystem 

k  =  the  number  of  components 

p.  =  the  reliability  of  the  i-th  component 

M.  =  the  exponent  of  the  i-th  component 

The  computer  program  modeled  a  system  that  had  8  sub- 
systems and  14  components.   System  reliability  (RS)  was 
determined  for  each  case  and  a  lower  confidence  bound  for  a   =  . 1 
and  a  =  .2  was  computed.   Random  numbers  were  drawn  from  a 
shuffled  random  number  generator  [Ref.  3].   Inverse  chi-square 
values  were  determined  using  the  international  mathematical 
and  statistical  library  (IMSL)  routine  called  MDCHI.   All 
computations  were  done  in  single  precision  arithmetic,  coded 
in  FORTRAN,  using  an  IBM  3  60  computer. 

A.   ZERO  FAILURE  VS  ONE  FAILURE  CASE 

An  examination  of  two  cases  revealed  a  shortcoming  and 
a  motivation  for  evaluating  the  modified  log-gamma  procedure. 
These  two  examples  are  considered  below. 
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Example  1. 

Let  k,  the  number  of  components  in  the  system,  be  14 
and  let  R.,  the  component  reliabilities,  all  equal  .99. 
The  sample  sizes  (mission  trials)  and  failures  for  each 
component  are  listed  in  Table  I.   The  lower  90%  confidence 
limit  on  system  reliability  is  desired. 

Table  I 

Component 
1   2   3   4   5   6   7   8   9   10  11  12  13  14 
M. :   #  mission  trials   10  10  10  10  500  10  10  10  10  500  10  10  10  10 

f.:      #  failures  00000000000000 


When  the  Zf .  =  0  the  modified  log-gamma  procedure  defines 


N*,  the  effective  sample  size,  as 


N*   =    k  -  (3.2) 


i=l   l 


For  the  data  given  in  the  table  above  N*  is  equal  to  11.62  8. 
For  this  procedure  the  lower  1-a  confidence  bound  R(l-a)  is 
computed  according  to  a  binomial  confidence  bound  based  on 
zero  failures  out  of  N*  trials.   The  value  obtained  for 
R(l-a)  was  .820. 


16 


Example  2. 

Let  sample  sizes  and  f.  (failures  for  each  component) 
be  given  in  Table  II.   Again  the  lower  90%  confidence  limit 
on  system  reliability  is  desired. 

Table  II 

Component 
1   2   3   4   5   6   7   8   9   10  11  12  13  14 

M. :   #  mission  trials   10  10  10  10  500  10  10  10  10  500  10  10  10  10 

l 

f.:      #  failures  00001000000000 

When  Ef.  ^  0  the  modified  log-gamma  method  solves  for 
R(l-a)  (the  lower  confidence  bound)  by  the  following  procedure. 
Let 


fi 
R.   =   1  - 


i         n . 

l 


/ 


Rc   =   .998  ,  (3.3) 

3 


R.   =   1   for   i/   5 
l 


be  the  point  estimate  of  the  i-th  component  reliability. 
Then 


R   =    n   R.   =   .998  (3.4) 

i=l   1 
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is  the  point  estimate  for  system  reliability.   Define 


r   =   r1^   =   .99986  (3.5) 


V   =   (1-R)   J      —      =       .000172      (3.6) 
1=1   l 


where  V  estimates  the  variance  of  -In  R. 
Let 


"  2 
L*   =   ("1^  R)    =   .02328  (3.7) 

V 


and  define 


L   =   L*  +  2.25   =   2.27328         (3.8) 


where  2.25  is  the  correction  term  and  L*  is  the  method-of- 
moments  estimate  of  the  shape  parameter.   Then  the  lower 
1-a  confidence  bound,  R(l-a)  is  computed  by  solving 

,(2L/x.2L  a} 
R(l-a)   =   R      /lj,a  (3.9) 

2 
where  Xo-P    is  t^le  l°wer  a  quantity  of  the  chi-square 

distribution  with  2L  degrees  of  freedom.   In  example  2  R(l-ot) 

is  equal  to  .99  3 

These  two  examples  have  shown  the  shortcoming  of  this 

method.   The  lower  confidence  bound  for  one  failures  is  higher 

than  the  lower  confidence  bound  for  zero  failures. 
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B.   SIMULATION  RESULTS 

The  lower  confidence  bound  values  obtained  for  the  twelve 
cases  studied  have  been  listed  in  Table  III.   RS  is  the  system 
reliability,  ACV  is  the  actual  confidence  value  computed  by 
the  modified  log-gamma  method  and  R(l-a)*500  is  the  percentile 
value  of  the  500  ordered  R(l-a)  estimates  for  a  =  .1  and 
a  =  .2.   N(I),  RI(I)  and  W(I)  are  the  respective  sample  sizes, 
reliabilities  and  weights  assigned  to  each  case. 

For  example,  in  case  number  3  the  number  of  components 
k,  is  equal  to  14  with  the  sample  sizes  equal  to  50   -Tor 
i  t^  5  or  10  and  250  for  i  =  5  or  10.   The  reliabilities  of 
each  component  is  .99  and  the  8  weights  are  all  equal  to 
.125.   System  reliability,  R_ ,  was  computed  to  be  .816  and 
for  a   =    .1  the  450-th  value  in  the  ordered  500  LCL  estimates 
was  .895.   The  Rg  value  of  .816  was  the  35th  of  the  500 
ordered  LCL  estimates  yielding  an  actual  confidence  level  of 
7.8%.   Likewise  for  a  =  .2  the  400-th  value  in  the  ordered 
500  LCL  estimates  was  .898.   The  Rg  value  of  .816  was  the 
13-th  of  the  500  ordered  LCL  estimates  yielding  an  actual 
confidence  level  of  2.8%.   In  only  one  case  (case  8)  did  the 
actual  confidence  value  approach  that  of  the  system  reliability 
as  a  lower  bound. 

An  examination  of  the  MLG  (modified  log-gamma)  procedure 
questioned  the  inclusion  of  the  correction  term  2.25.   Addi- 
tional simulations  were  run  on  the  same  twelve  cases  when  this 
correction  term  was  removed  and  the  degrees  of  freedom  bounded 
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below  by  1.0.   The  results  obtained  from  this  modification, 
while  an  improvement,  were  still  far  from  providing  accurate 
lower  bounds  on  the  system.   The  values  determined  from  these 
runsaJre  listed  in  Table  IV.   ACV  values  of  100%  indicate  that 
the  system  reliability  was  greater  than  all  500  estimates. 

It  would  appear  that  in  order  to  generate  more  estimates 

2 
less  than  RS  the  exponent,  2L/y  ~   /  needs  to  take  on  larger 

r  A2L,ct 

values.   Adding  a  constant  term  such  as  2.25  yields  more  values 
for  R(l-a)  that  are  greater  than  RS .   Indeed,  Tables  III  and 
IV  did  show  this  to  be  the  case.   As  the  exponent  becomes 
larger  (the  chi-squared  value  smaller)  the  confidence  level 
decreases.   The  estimate  for  L  used  in  generating  the  values 
listed  in  Table  IV  seem  more  accurate  when  used  in  the 
modified  log-gamma  procedure. 

This  modification  still  left  much  room  for  improvement. 
A  closer  reivew  of  the  MLG  method  pointed  to  the  estimate  of 
the  shape  parameter  as  a  possible  cause  of  the  extreme  results. 
Since  Z  =  -In  R  its  distribution  was  approximated  by  a  two- 
parameter  gamma  distribution.   Then 


L-l   Lz/ln  k 

f(z)   =   - ■  e      T ,   z  >0       (3.10) 

(-lnR_)  L  r(L) 


In  R 
where  L  and  (-  — - — )  are  the  parameters.   Then 


E(z)   =   L  •  ("1"  R)   =   -  In  R  (3.11) 
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and 


2 
Var(z)   =   M"1!?  R)2   =   ±^3  (3.12) 


Note: 


l   =  %g*  =   IgMl!  ,3.13) 

Var(z)      Var(z) 


The  proposed  estimator  L  for  L  is 


Z2 

L    =  -^—  (3.14) 

Var(z) 

2  ^  2 

fE  (  z )  1  ^    [E(z)l 

Since  L  =  ^ — ^rM —   it  would  appear  that  L  =  - — ysr— - —  would 
Var(z)  ^  „A,  x 

Var(z) 

be  a  better  estimator  for  L.   Since 


[E(z)  ]2   =   E(z2)  -  Var(z) 


we  have 


^       A 

?      E(z  )  -  Var(z)  ,-  ,_ 

L   ~   VaHzl (3'15) 


2  2 

and  since  z   is  unbiased  for  E(z  )  we  get 


t      =   z2  "  Var(z)   =   _J^ !         (3.16) 

Var  ( z )         Var I z ) 
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Note  that  this  is  a  departure  from  L,  in  the  proposed 
method.   Thus  the  shape  parameter  L  can  be  estimated  by 
Eq.  3.16  above.   This  estimate  is  different  from  the 
original  version  of  the  MLG  method. 

Substituting  this  new  value  for  L  and  bounding  the  degrees 
of  freedom  by  1.0,  so  as  not  to  obtain  a  negative  value,  the 
results  show  a  little  more  improvement.   The  results 
obtained  from  this  second  modification  are  listed  in 
Table  V. 
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IV.   CONCLUSIONS 

Additional  simulations  on  many  more  cases  would  be 
required  to  determine  the  particular  conditions  under  which 
this  modified  log-gamma  method  is  reasonably  accurate.   For 
the  cases  examined  here  the  proposed  procedure  remains 
suspect  in  estimating  lower  confidence  bounds  on  system 
reliability. 
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APPENDIX 


AA  CORRECT  ICIV    TERM    EQUAL    TO    2.25     IN    THE    MODIFIED    LOG- 

GAMMA    METHOD 

A8  VARIABLE    THAT    STORES    THE    CIFFERENCE    BETWEEN    PS 

(SYSTEM    RELIABILITY)     AND    RR(^OO) —    THE    80"TH 
PERCENTILE    POINT    WHEN    ALPHA  =  0.2 

ABS  ABSOLUTE    VALUE 

AC  VARIABLE    THAT    STORES    THE    CIFFERENCE    BETWEEN    ^S 

(SYSTEM    RELIABILITY)     AND    R(450) —    THE    90TH 
PERCENTILE   POINT    WHEN    ALPHA  =  0.1 

ALOG  NATURAL    LOGARITHM    SLBROUTINE 

ALPHA  VARIABLE    ASSIGNED    A    VALUE    OF    0,1 

ALPHAA  VARIABLE    ASSIGNED    A    VALUE    OF    0.2 

AM  ARRAY    THAT    STORES    THE    EXPONENTS M    SUB    I 

8LHAT  VARIABLE    THAT    STCRES    THE    L    HAT    VALUE 

CA  ACTUAL    CCr>FIDENCE    LEVEL     FOR    ALPHA=0.1 

CALL  FORTRAN    CODE    FOR    ACCESSING    SUBROUTINES 

CB  ACTUAL    CONFIDENCE    LEVEL    FOR    ALPHA=0.2 

CONTINUE  FORTRAN    CCDE    TO    CLOSE    EACH    00    LCCP 

CA  DUMMY    VARIABLE    USED   TO    DETERMINE    THE    ACTUAL 

CONFIDENCE    LEVEL 

DDF  CEGPEES    OF    FREEDOM 

DIMENSION  FORTRAN    CCDE    REQUIRED    FOR    DIMENSIONING    ARRAYS 

DC  FORTRAN    CCDE   US  EC   TO    BEGIN    LOOPS 

DUM  CUMMy    VAR  IABLE 

EA  DUMMY    VARIABLE    LSED    TO    DETERMINE    ACTUAL 

CONFIDENCE   VALUE 

EFFN  EFFECTIVE    SAMPLE    SIZE 

END  FORTRAN    CCDE    REQUIRED  TD     END    PROGRAM 

EXP  EXPONENTIAL    SUBROUTINE 

FA  DUMMY    VARIABLE     LSED    TO    DETERMINE    THE    ACTUAL 

CONFICENCE   VALUE 

FLOAT  FORTRAN    CCDE    USED    TO    CHANGE    INTEGERS    to    CECIMAL 

VALUES 

FORMAT  FORTRAN    STATEMENT 

GA  DUMMY    VARA 

GA  DUMMY    VARIABLE     USED    TO    DETERMINE    THE    ACTUAL 

CONFI  CENCE   VALUE 

GO  FORTRAN    CODE    USED    IN    THE    — GO    TO — STATEMENT 

HISTG  SUBROLTINE    WHICH    GENERATES    A    HISTOGRAM    CF    THE 

DATA     AND     SUMMARY    STATISTICS 
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INDEX    VAR  IA3  LE 

ERROR    VARIABLE     IN    SUBROUTINE    MCOI 

VARIABLE    USED    TO    STORE    THE    NUMBER    IF    FAILURES 
(ALSO    PART    OF    THE    FORTRAN       — IF —   STATEMENT) 

INDEX    VAR  IABLE 

VARIABLE     THAT    STORES    THE    INITIAL    VALUE     FCF    CALLING 
RANDOM    NUMBERS 

INDEX    VARIABLE 

VARIABLE    THAT    STORES    THE    NUMBER    OF    FAILURES 
PER    COMPONENT 

INC  EX    VAR  IABLE 

INDEX    VAR  IA3  LE 

VARIABLE    THAT    STCRES    THE    NUMBER    CF    COMPONENTS 
(ALSO    USED    AS    AN    INDEXING    VARIABLE) 

INDEX    VARIABLE 

VARIABLE    THAT    STORES    THE    NUMBER    OF    SUBSYSTEMS 
(ALSO   USEC    AS     AN    INDEXING   VARIABLE) 

COUNTER    VARIABLE 

INVERSE    CHI    SQUARE    SUBROUTINE 

INDEX    VAR  IABLE 

ARRAY    THAT    STORES    THE    K    SAMPLE    SIZES 

COUNTER    VARIABLE 

VARIABLE    THAT    STORES    THE    NUMBER    OF    CASES 

INDEX    VARIABLE 

SUBROUTINE    REQUIRED    FOR    RANDOM    NUMBER    GENERATION 

ARRAY    THAT   STORES    THE    UNIFORM    RANDOM    NUMBERS 

VARIABLE    THAT    STCRES    THE     4L°H/i    V  4  LU  E    CF     .1 

VARIABLE    THAT    STCRES    A    PO  IN""-     ESTIMATE 

VARIABLE    THAT    STORES    A    POINT    ESTIMATE 

ARRAY    THAT   STORES    THE    LOWER    CONFIDENCE    BOUNC 
VALUE    WHEN    ALPHA=0.1 

VARIABLE    THAT    STORES    THE    R(450)     VALUE 

VARIABLE    THAT    STCRES    R    BAR 

PvUMBER    CF    REENTRY     ECDI ES    °  ER    MISSILE 

CUMMY    VARIABLE    USED   TO    COMPUTE    RBAR 

FORTRAN     STATEMENT 

VARIABLE    THAT    STCRES    THE    INVERSE    3  F    THE    EFFECTIVE 
SAMPLE    SIZE 

RHAT  VARIABLE    THAT    STCRES    THE    SUM    OF    THE    WEIGHTED 

SUBGROUP    RELIABILITY    ESTIMATES 


37 


1 

IER 

If 

11 

I  SEED 

J 

JF 

JJ 

JM 

K 

KJ 

L 

MC 

MOCHI 

MM 

N 

NC 

NCASE 

NN 

CVFLOW 

P 

PD 

PN 

FR 

R 

RB 

RBAR 

RBN 

RCEN 

READ 

REFFN 

RI  ARRAY    THAT    STORES    THE    IMPUTED    RELIABILITY    VALUES 

RIHAT  ARRAY    T  HAT    STORES    THE    COMPUTED    RELIABILITY    VALUES 

RWEAN  VARIABLE    THAT    STORES    THE    MEAN    OF    THE     R    ARRAY 

RNUM  DUMMY    VARIABLE    USED    TO    COMPUTE     RBAR 

RR  ARRAY    THAT    STORES    THE    LOWER    CCNFIDENCE    EOUNCS 

WHEN    ALPHA=0.2 

RR8  VARIABLE  THAT    STCRES    THE    RR(400)    VALUE 

RRMEAN  VARIAELE  THAT    STORES    THE    MEAN    OF    THE    RR    AFRAY 

RRVAR  VARIABLE  THAT    STORES    THE    VARIANCE    OF    THE     RR    ARRAY 

RS  VARIABLE  THAT    STORES    THE    TOTAL    SYSTEM    RELIAEILITY 

RUHAT  ARRAY    THAT    STORES    THE    SUBGRCUP     RELIABILITY 

ESTIMATES 

RVAR  VARIABLE    THAT    STCRES    THE    VARIANCE    OF    THE     F    ARRAY 

S  ARRAY    THAT    STHRES    THE    VAR/COV    MATRIX 

SCR  VARIABLE    THAT    STORES    THE    STANDARD    DEVIATION    OF 

THE    R    ARRAY 

SDRR  VARIABLE    THAT    STORES    THE    STANDARD    DEVIATION    OF 

THE    RR    ARRAY 

SGRT  SUBROUTINE   THAT    SOLVES    SQUARE    ROOTS 

SRAND  SUBROUTINE    THAT    IS    THE    SHUFFLED    RANDOM    NUMBER 

GENERAT  OR 

STOP  FORTRAN    REQUIRED    CODE 

SUM  DUMMY    VARIABLE    USED    THROUGHOUT    ThE    PROGRAM 

TG  PART    OF   TFE    FORTRAN    — GO    TO —    STATEMENT 

V  ARRAY    THAT    STORES    4    VARIANCE    ESTIMATES    FOR 

COMPONENTS    11    THROUGH    14 

VFAT  VARIABLE     THAT    STORES    THE    VARIANCE    ESTIMATE    FCR 

-LN(RFAT  ) 

VN  VARIABLE    THAT    STORES    THE    VARIANCE    ESTIMATE    FOR 

-LN(PN) 

VR  VARIABLE     THAT    STORES    THE    VARIANCE    ESTIMATE    FOR 

-LN(PP) 

VX  DUMMY    VARIABLE     USED    IN    THE    MOCHI     SUBROUTINE 

VY  DUMMY    VARIABLE    USED    IN    THE    MCCHI    SUBROUTINE 

W  ARRAY    THAT    STORES    THE   WEIGHT  EC    VALUES     OF    EACH 

SUBSYSTEM 

WRITE  FCRTRAN    STATEMENT 

Z  DUMMY    VARIABLE    USED    TO    DETERMINE    RBAR 
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APPENDIX     B 


DIMENSION    R(  5C0)  ,RI  (5  0)  ,W(5  0)  |N(50J  ,PIHAT  (5  0)  VA»  (50  ), 

1PUHAT  (50),V  (4),S(  8,  8)  ,RR(  500)  ,P(  5CC) 
CALL    CVFLGW 
NCASE    =    0 
ISEED    =     134869 

C 

c 

C  PEACING     IN    THE    INPUT    PARAMETERS    K  ,L  ,RBN,A  ,  A  LPHA    AND 

C  ALPHAA 

C 

c 

READ    (5,330)     K,  L,  RBIS,  A  A,  ALPHA  ,  AL  P  FAA 
C 
C 

C  PEACING    IN       N(I) THE   NUMBER    OF    R/NCQM    nhjmeERS    P  ER 

C  COMPONENT 

C 

c 

20    REAC    (5,340)     (N(I),I  =  1,K) 
C 
C 

C  READING    IN    THE   COMPONENT/FUNCTION    RELIABILITIES 

C 
C 

REAC    (5,350)    (RI(  I  ),  1=1, K) 

NCASE    =    NCASE+1 

IF     (  ISEED. GT.  134869)    GO    TO    30 
C 
C 

C  READING    IN    THE    EXPONENTS   M    SU8    I 

C 
C 

REAC    (5,360)     (AM(  I),  1=1, K) 
C 
C 

C  RECCING    IN    THE   WEIGHTS     FOR    EACH    SUB  SYSTEM /G  RO  LP 

C 
C 

30     PEAC     (5,370)     (W(  I),  1=1, L) 
C 
C 

C  STARTING    THE    NAIN    LCCP    FOR    500    SIMULATIONS 

C 
C 

c 

DO     190  1  =  1,5GC 
C 

C  LOOPING    FOR    EACH    COMPONENT 

C  AND    DRAWING    THE    RANDOM    NUMBERS 

C 

IF     =    C 
C 

OG    50    J=l  ,K 

JJ    =    N(  J) 

C4LL    SRANC    (  ISEED, o,JJ  ) 

JF    =    0 
C 

CO    40    JM=1,  JJ 

IF     <P(JM)  .GT.RK  J)  )     JF  =  JF  +  1 
40    CONTINUE 


IF    =    IF+JF 
C 
C  VARIABLE    "IF"    COUNTS    THE    FAILURtS 


C 


RIH/TU)    =     l.-(FLOAT(  JFJ/FLOATl  JJ)  ) 

50    CONTINUE 
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c 

c 


JF     (IF.E0.0  )    GC   TO    170 
C 

C  COMPUTING    THE    POINT    ESTIMATES 

C 

PN    =    1. 

PR    =    1. 

CO    60    J=lt  5 
FR    =     (RIHAT(J)**AM(J) )*PR 
PN    =    (P.IHATC  J+5)**AM<  j+5)  )*PN 
60     CONTINUE 
L 

c 
c 

C  COMPUTING    THE    SUBGRCUP    RELIABILITY     ESTIMATES 

C 

C 

RUHAT(l)     =   PR*RI(11J*RI(12J 

RUH6T(2)    =    PN*RI  (  11)*RI  (12) 

RUH£T<3)     =    PR*Ri(13  J*RIU4) 

RUHAT(4)    =    PN*RI<13)*RI(14i 

RUH*T(5J    =    PR*R  I(  13  J*o  1(12) 

RUHATC6)     =    PN*RI (13  )*RI (12) 

RUHAT(7)     =    PR*RI(  11)*RI  (  14) 

PUhAT(8)     =    PN*RI( 11  )*RI(  14) 
C 

C  COMPUTING    RHAT    AND    CALLING    IT    BY    THE    SAME    NAME RHAT 

C 

RHAT    =    0. 
C 

OC    7C    J=l tL 

RHAT    =    ( W( J)*RUHAT( J) J+RHAT 
70     CONTINUE 
C 
C 
C 

C  ESTIMATING    THE    VARIANCE    OF    -LMRH4T) VF4T 

C 
C 
C  STEP     l:       DETERMINING    RBAR 


C 


C 
C 


RNUM    =    0. 
RDEN    =    0. 

CO    80    J=1,K 

RNUN    =     (-4L0G(RIHAT (J) )*AM(J ) ) +RNUM 

RDEN    =    AM( J )+RDEN 
80     CONTINUE 
C 

Z   =     (-RNUM) /RDEN 

IF     (Z.LT.O)    GC   TO   90 

RBAR    =   EXPC  Z) 

GO    TO    100 
90     RBAR    =    l./EXP(ABS ( Z) ) 
C 
C  STEP    2:        DETERMINING   THE    VARIANCE    ESTIMATES 


r 


C 

c 


100    VR    =    C. 
VN    =   0  . 

CO    110    J=l,  5 

VR    =     (AM( J)**2/FL0AT(N(J ) ))+VR 
VN    =    (  AM(  J+5)  **2/FL0AT(N<  J+5)  )  )  +  VN 
110     CONTINUE 

VR    =    (  l.-RBAR  )*VR 
VN    =     ll.-RBAP)*VN 

DO     120    J=lr4 

V(J)     =    (1  .-RB4R)/FLCAT  (N(J+10  )  ) 
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c 
c 
c 
c 
c 
c 


c 
c 

c 
c 


c 
c 
c 
c 

c 
c 


c 

c 
c 


120    CONTINUE 


STEP  3 


FINAL  SOLUTIONS  FOR  COVARIANCE  ESTIMATFS 


S(L 
S(2 
S(3 
S(  4 
S(5 
S(6 
S(  7 
S  (8 
S(l 
S(  1 
S(l 
S(  1 


,11 
,2) 
i3) 
,4) 
,5) 
,6) 
,7) 
,8} 
,2) 
,3) 
,4) 
,5) 


S( 1,6)    = 


S(l 

Si  1 

S(2 

S(2 

S(  2 

S(2 

S(2 

S(2 

S(3 

S(  3 

S(3 

S<3,7) 

S(  2,  6) 

S(4,5) 

S(4,6) 

S(4,7) 

S<4  ,8  ) 

S(  5 

S(  5 

S(5 

S<6 

S  (6 


,7) 
,8) 
,3) 
,4) 
,5) 
,6) 
,7) 
,6) 
,4) 
,5) 
,6) 


,6) 
,7) 
,8) 
,7) 

,8J 


S(7,8)    = 


VR  +  V 
VN+V 
VR+V 
VN+V 
VR+V 
VN+V 
VR  +  V 
VN+V 

vm 

VR 

0. 

VR+V 

V(  2) 

VR+V 

V(  1) 

0. 

VN 

V(  2) 

VN+V 

V(  1  J 

VN+V 

V(3  J 

VR+V 

V(3> 

VR+V 

V(  4) 

V(3  J 

VN  +  V 

V(  4) 

VN+V 

V(3) 

VR 

0. 

0. 

VN 

V(l) 


<1)+V(2) 
(  1)  +  V(  2) 
(3)+V  (4) 
<3)  +  V(4) 
(3)+V<  2) 
(3  1+V  (2  J 
(  1)  +  V(  4) 
(1  J+V(4) 
♦  V(2  J 


(2) 
(1) 


(2) 

(  1) 

•»V(4  ) 
(3) 

(4) 


(3) 

(4) 

+  V(2) 


+V(4) 


FILLING    IN   XHE    REST    OF    THE    VAR/COVAR    MATRIX 


DO    140    MM=1,L 


CD    130    NM=1,L 
SINN, MM)     =    S(MM,NN) 
130    CONTINUE 

140    CONTINUE 


SOLVING   THE    OVERALL    EQUATION    FOR    VHAT 

VHAT    =   0 . 

CO     160    J=l,  L 

DO     150    KJ=1,L 

VHAT    =    W(  J)*W(KJ)*«*UhAT(J)*RUHAT{  KJJ*S<  J,KJ)  +  \HAT 

150    CONTINUE 
160     CONTINUE 

VHAT    =    VHAT/(  RHAT**2) 
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c 

c 
c 
c 

c 
c 
c 
c 
c 
c 
c 


c 
c 
c 
c 

c 
c 


c 
c 

r 


c 
c 
c 
c 
c 


c 

c 


c 
c 

c 
c 
c 
c 


COMPUTING    LFAT 


eLF*T=((-ALOG(RHAT)  )**2/VHAT)-  l.Q 


COMPUTING    ThE    DEGREES    OF    FREEDOM-DDF    AMD    SOLVING    FOR 
OF     (1-ALPHA)     VHEN    THE    SUM    OF    THE    FAILURES     CCES    NnT 
EQUAL    ZERO 


DDF    =    2.*BLHAT 

IF    (CCF.LT.l.O)    DOF-1.0 

PD    =    .  1 

CALL    MDCHI    (PCtODFt  VX,  IER  ) 

PD    =    .2 

CALL    MDCHI    ( P D,DDF , VY ,1 ER ) 

R(I J    =    R  FAT** ( (DDF  J/VX  ) 

RR(I)    =    RHAT**((DDF)/VY ) 

GO    TO    190 


COMPUTING  RELIABILITY  ESTIMATES  WHE>  THE  SUM  CF  THE 
FAILURES  IS  GREATER  THAN  ZERO 


170    SUM    =    0. 

DO    130    II =1  ,K 

SUM    =    SUM-M  l./FLOAT(N( II) J ) 
180    CONTINUE 

EFFN    =    FLOATCKJ  /SUM 

REFFN   =     1  ./EFFN 

R(I  )     =    ALFHA**REFFN 

RR(  I)    =    ALPHAA**REFFN 

GO    TO    190 
190    CONTINUE 


CALL    HISTG    (R ,500»  0) 
CALL    HISTG     (RF,500t  0) 


COMPUTING    THE    TOTAL    SYSTEM    RELIABILITY 


DUM    =    0. 

DC    210    J=1,L 
SUM    =    1. 

DO    200    1=1, K 
SUM    =    SUM*(RI(  I)**AM(  I  )  ) 
200     CONTINUE 

DUM    =    (W( J) *SLM)+DUM 
210     CONTINUE 

PS    =    DUM 


AC    =    PS -R (450) 
A3     =    RS-RR(400) 
R=*    =    R( 450) 
PRB    =    RR(400) 


COMPUTING   ThE   SAMPLE    VARIANCE    F3R     ^LPHA=.l — SCR 
AND    ALPHA=.2 SDRR 
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220 


230 


SUM    =    0. 

CU  M    =    0  . 

01    220    1=1,  500 
SUM    =    R(I J+SU f 
OUM    =    RRUJ+DLM 
CONTINUE 

RMEAN    =    SUM/5C0. 
RRMEAN   =     CUM/ 500. 
SUM    =    0. 

OJM    =    0. 

03    230    1=1,500 

SUM    =    (R(  I  )-RMEAN)**2+SUM 

£'JP    =    (RR(IJ-PRMEAN  )**2+DUM 

CONTINUE 

PVAF    =    SUM/499. 
RRVAR    =    DUM/4<5. 
SDR    =    SQRT(RVAR) 
SDRP    =   SQRT  (R.RVAR) 


C 
C 

c 
c 


c 
c 


240 
250 

260 
270 


DETERMINING    TFE    ACTUAL    CONFIDENCE    V4LUE 


MC    =    0 

CC    250    1=1,500 

IF     (RS.GT.Rd  ))    GO    TC    240 

DA    =    R(I )-RS 

EA    =    RS-R(MC) 

IF     (DA.LT.EA)     GO    TO    260 

CA    =     FLOAT(MC  J/500  . 

GO    TO    270 

MC    =    MC+1 

CONTINUE 

CA    =    FLOAT  (  D/500. 
CA    =    CA*100. 


NC    =    0 


C 

c 
c 
c 

r 


280 
250 

300 
310 


DO     290    I =1,5CC 
IF     (RS.GT  .RR(  I))     GO    TO     280 
FA    =    RR(I)-RS 
GA    =    RS-SR(NC) 
IF     {FA.LT  .GA)    GO    TO     300 
C3    =    FLOAT!  NC  )/5  0O. 
GO     TO    310 
NC    =    NC+1 
CONTINUE 

C8    =    FLCAT(IJ/500. 
C3    =    CB*1C0. 


PRINTING    THE    FINAL    RESULTS    FOR    EACH    CASE 


WRITE 

V»RITE 
WI  TE 
W  R  IT  E 
ViRI  TE 
ISEEO 
IF 
GO 


NCASE 


310     STOP 


(6,380) 

(6,390) 

(  6,  400) 

(6,410  J 

(  6  ,  42  0  ) 

=     ISEED  +  27 
(NCASE. EC. 8)    GO 
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